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Abstract. - The spin-polarized transport through two- level quantum dots weakly coupled to 
ferromagnetic leads is considered theoretically in the Coulomb blockade regime. It is assumed 
that the dot is doubly occupied, so that the current flows due to cotunneling through singlet 
and triplet states of the dot. It is shown that transport characteristics strongly depend on 
the ground state of quantum dot. If the ground state is a singlet, differential conductance 
(G) displays a broad minimum at low bias voltage, while tunnel magnetoresistance (TMR) is 
given by the Julliere value. If triplet is the ground state of the system, there is a maximum 
in differential conductance at zero bias when the leads are magnetized in antiparallel. The 
maximum is accompanied by a minimum in TMR. The different behavior of G and TMR may 
thus help to determine the ground state of the dot and the energy difference between the singlet 
and triplet states. 



Introduction. Transport properties of quantum dots coupled to ferromagnetic leads 
have been a subject of thorough studies since a few years [1,2]. The considerations concerned 
both the strong coupling regime where the Kondo physics emerges [3], as well as the weak 
coupling regime. In the weak coupling regime most of theoretical works addressed the problem 
of spin-dependent sequential transport through quantum dots hosting single orbital level [4] . 
In the sequential tunneling regime, if the applied bias voltage exceeds a certain threshold 
voltage, electrons tunnel one by one through the system, otherwise the current is suppressed 
leading to the Coulomb blockade effect. Although in the Coulomb blockade regime the se- 
quential transport is suppressed, the current can be still mediated by higher-order tunneling 
processes such as cotunneling which involves correlated tunneling of two electrons via virtual 
states of the dot [5] . Spin-polarized cotunneling transport has been addressed very recently. 
It has been shown that when the dot is singly occupied, a zero-bias anomaly appears in dif- 
ferential conductance when the magnetizations of the leads are aligned in antiparallel [6,7]. 
These considerations were performed only for single-level quantum dots. In real systems, how- 
ever, usually more than one energy level participate in transport, which can lead to further 
interesting effects [8]. 

The goal of this paper is to address the problem of spin-dependent transport through two- 
level quantum dots coupled to ferromagnetic leads in the cotunneling regime. More specifically, 
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we analyze the case when the dot is doubly occupied at equilibrium and the system is in the 
Coulomb blockade regime. In this case, depending on the value of the exchange interaction, 
either the singlet or triplet state will be favored. We show that the differential conductance 
(G) and tunnel magnetoresistance (TMR) display a strong dependence on the ground state of 
the system. If the ground state is a singlet, G exhibits a broad minimum at low bias voltage 
in both magnetic configurations of the system, while TMR takes the Julliere value [9] . On the 
other hand, if the ground state is a triplet, there is a maximum in differential conductance at 
zero bias in the antiparallel configuration, accompanied by a minimum in TMR. The different 
behavior of both differential conductance and TMR provides thus a handle to determine the 
ground state of the system. 

The systems discussed in this paper may be realized experimentally for example in metallic 
single- wall carbon nanotube (SWCNT) quantum dots contacted to ferromagnetic leads [10]. 
As shown recently, the shell filling of SWCNT quantum dots exhibits four-electron periodicity 
when sweeping the gate voltage [11]. By changing the gate, one can tune the system to the 
Coulomb blockade regime where the dot is doubly occupied at equilibrium. 

Model and method. - We consider spin-polarized cotunneling transport through quantum 
dots with two orbital levels coupled to ferromagnetic leads. It is assumed that the magnetic 
moments of the leads can form either parallel or antiparallel magnetic configuration, see fig-H 
The Hamiltonian H of the system involves four terms, H = Hl+Hr+Hd+Ht- The first two 
terms describe the noninteracting itinerant electrons in the leads, H r = ^ fc(T £ r kac\, k(J c r k a , for 
r = L,R, corresponding to the left and right leads, with £ rk cr being the energy of an electron 
with wave number k and spin a in the lead r, and c\ ka (c rka ) denoting the respective creation 
(annihilation) operator. The quantum dot is described by the Hamiltonian [12] 

H D = y^jtijc + U/.n^riji + U'^n la n 2 a' - J A.a d H3 d \y d i5^ a ^s , (1) 

j(7 j go' a fly 6 

where nj a — dj^dja and dt {dj„) is the creation (annihilation) operator of an electron with 
spin a on the jth level (j — 1,2), while Ej is the corresponding energy. On-level and inter- 
level Coulomb interaction between electrons is described by U and U', respectively. The last 
term in Hp corresponds to the exchange energy due to the Hund's rule, with J being the 
respective exchange coupling and a denoting a vector of Pauli spin matrices. To simplify 
further discussion of numerical results we assume U' — U. For the energy of the first orbital 
level we write E\ = e and e 2 — £ + Se, where Se is the splitting of the dot orbital levels. 
For example in metallic SWCNT quantum dots, 5e would correspond to the energy mismatch 
between the two subbands of the nanotube [11]. The mismatch can vary with the boundary 
conditions at the nanotube-electrode interface. 

The last term of H describes the tunneling processes between the dot and electrodes 

Ht= J2 E (TrdkJia + T^d)^) , (2) 

r— L,R jka 




Fig. 1 - Schematic of a quantum dot coupled to ferromagnetic leads. The magnetizations of the leads 
can form either parallel or antiparallel magnetic configuration. The system is symmetrically biased. 
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with T r j denoting the tunnel matrix elements between the lead r and jth level. Coupling 
of the jth dot level to external leads can be described by T 



with p a r being 

the spin-dependent density of states in the lead r. By introducing the spin polarization of 
lead r, p r = (p+ — p~)/(p+ + Pr)> the couplings can be expressed as, T^f = T r j(l ±p r ), 
with T r j = (r+ + r ? T J )/2. Here, and T~ describe the coupling of the jth level to spin- 
majority and spin-minority electron bands, respectively. We assume that the system is coupled 
symmetrically to the leads T r j =T/2 and j?l = Pr = P- As reported in ref. [13], typical values 
of T in the weak coupling regime are of the order of tens of peV. 

In order to find the current flowing through the dot in the Coulomb blockade, one has to 
calculate the respective cotunneling rates. Within the second-order perturbation theory, the 
rate of a cotunneling process from lead r to lead r' which changes the dot state from \x) into 



|x') is given by 



2tt 

T 



E 



m\H T \$ v )($ v \H T \$x) 



5{e l - s f ), 



(3) 



with £j and Sj denoting the energies of initial and final states, |$*) being the state of the 
system with an electron in the lead r and the dot in state |x) 5 whereas is a virtual state 
with e v denoting the corresponding energy. Among different cotunneling processes one can 
distinguish the single-barrier (r = r') and double-barrier (r / /) cotunneling as well as spin- 
nip (x x') and non-spin-flip (x = x') cotunneling. The spin-flip processes change the state 
of the dot, whereas the non-spin-flip processes are fully coherent, i.e. they do not change the 
dot state. The current flows through the system due to double-barrier cotunneling processes. 
On the other hand, the single-barrier processes do not contribute directly to electric current, 
however, they can change the dot occupations, and this way also the current. The details of 
how to determine different cotunneling rates can be found in ref. [5]. 

The cotunneling current flowing through the system from the left to right lead is given by 



XX' 



(4) 



where P x denotes the corresponding occupation probability. 



The probabilities ± x 



-7 X 7 

I rr 



• 7 

I r 



1 v 



P x can be 
together with the 



found from the master equation, = Ylrr' 
normalization condition P x = 1. 

We consider the case when the dot is doubly occupied at equilibrium and the system is 
in the Coulomb blockade regime. There are then six different two-particle states in the dot 
possible, these are three singlets \S = 0,M = 0)i = (| T)| 4) - I 1)1 T))/V2, |0,0) 2 = | U)|0), 
|0,0> 3 = 1 0> I U>, and three triplets |1,0) = (| T)| l> + I 4)1 T»/A |1,1) = I T)| T> and 
1 1,-1) = | |)| |), where the first (second) ket corresponds to the first (second) orbital level 
of the dot. In the case of finite level spacing, 6e > fceT, T, and J < Se, the lowest singlet 
state is |0, 0)2- In the following, we show that transport characteristics strongly depend on 
the ground state of the system. It is therefore useful to introduce the difference between the 
energy of the lowest lying singlet (es) and triplet (et) states, Ast = £s — £t = J — 5e. 

Results and discussion. - Figure |5] presents the bias voltage dependence of differential 
conductance G in the parallel and antiparallel magnetic configuration for several values of 
Ast = J — Se. First of all, one can see that the transport characteristics show a distinctively 
different behavior depending on Ast, i-e- the ground state of the quantum dot. For Ast < 0, 
the ground state of the dot is a singlet, |0, 0) 2 , whereas for Ast > 0, the ground state is a 
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Fig. 2 - (color online) The differential conductance in the parallel (solid line) and antiparallel (dashed 
line) magnetic configuration of the system as a function of the bias voltage for several values of 
A ST = J - 5e as indicated: (a) A ST = -5I\ (b) A ST = -2.5I\ (c) A ST = 0, (d) A ST = 2.5r. The 
other parameters are: k B T = 0.5r, e = -60T, 5s = 5T,U = U' = 40r, and p = 0.5. 



triplet, which is three-fold degenerate, |1, 0), |1, 1), |1, —1). On the other hand, for Ast = 0, 
the dot is in a mixed state and the occupation of singlet and each triplet is equal at equilibrium 
and given by 1/4. 

First, let us discuss the case of Ast < shown in fig. [21 a) and (b). In this situation, at low 
bias voltage the dot is occupied by two electrons on the lowest energy level, |0, 0)2 = | tl)|0), 
and the ground state is singlet, 5 — 0. Because for Ast < and \eV\ < |Ast|> -P|o,o} 2 ~ 1> the 
current is mediated only by non-spin-flip cotunneling processes. However, once |eV| > |Ast|> 
the triplet states start participating in transport (in the case of Ast = — 5r also the state 
|0, 0)i contributes), leading to an increase in the differential conductance at \eV\ « |Ast|- 
In other words, the suppression of cotunneling through 5 = 1 states gives rise to a broad 
minimum in differential conductance which is present in both magnetic configurations of the 
system, see fig. |!Ja) and (b). The width of this minimum is determined by the splitting 
between the singlet and triplet states, 2|Ast|- By measuring the width, one can thus obtain 
information about the energy difference between the 5 = and 5=1 states. 

In the case when the current flows only due to non-spin-flip cotunneling through the 
state |0,0)2, one can find approximative formulas for the differential conductance. At low 
temperature and for Se ^> kBT, G in the parallel configuration can be expressed as 



(e + U) 2 {e + 5e + 2Uy 



(5) 



on the other hand, for the antiparallel configuration one finds, G^p° = (1 — p 2 )Gp =0 / (1 +p 2 )- 
These expressions approximate the minimum in differential conductance at low bias voltage 
shown in fig. Eta)- The analytical formulas for the general case when more states participate 
in transport are too cumbersome to be presented here. 
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Figure Pic) displays the differential conductance in the case when the singlet and triplet 
states are degenerate, Ast = 0. First of all, it can be seen that G exhibits a distinctively 
different behavior in both magnetic configurations. In the parallel configuration the differ- 
ential conductance has a smooth parabolic dependence on the bias voltage, whereas in the 
antiparallel configuration there is a maximum at the zero bias voltage. This effect bears a 
resemblance to that found in the case of singly occupied one-level quantum dots [6,7]. Here, 
however, the mechanism leading to the maximum is different - the zero-bias peak appears 
due to cotunneling through singlet and triplet states of the dot. In the case of Ast = and 
at low bias voltage, all the four dot states, i.e. |0, 0)2, 1 1, 0) , |1,— 1), |1, 1), participate in 
transport on an equal footing. Consequently, the current flows due to both spin-flip and non- 
spin-flip cotunneling processes. To understand the mechanism leading to the zero-bias peak, 
one should bear in mind that in the antiparallel configuration the spin-majority electrons of 
one lead tunnel to the spin-minority electron band of the other lead. For example, for positive 
bias voltage (electrons tunnel then from the right to left lead), the spin-| electrons can easily 
tunnel to the left lead (the spin-| electrons are the majority ones), while this is more difficult 
for the spin-| electrons (they tunnel to the minority electron band). Thus, with increasing 
the bias voltage, the occupation of state |1, —1) (| |)| J,)) is increased, while the occupation of 
state |1, 1) (| |)| t)) decreases. The unequal occupations of these triplet states lead thus to a 
nonequilibrium spin accumulation in the dot, P\\ y \) — P\\ —\) < 0, which is shown in fig. 01 It 
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Fig. 3 - (color online) Nonequilibrium spin accumulation in the dot, .FVi) — Pm-i), as a function of 
the bias voltage for Ast = —5, —2.5, 0, 2.5F. The other parameters are the same as in fig. [21 

is further interesting to note that in the antiparallel configuration the possible non-spin-flip 
cotunneling processes are proportional to T^T^ and r^T^, whereas the spin-flip cotunneling 
is proportional to T^T^ and r^T^. It is clear that the fastest cotunneling processes are the 
ones involving only the majority spins, i.e. r^r^. However, because of nonequilibrium spin 
accumulation, with increasing the bias (V > 0), the dot becomes dominantly occupied by 
majority electrons of the right lead, P\x,-i) > -P|x,i> 5 an d the processes proportional to r^T^ 
are suppressed. As a consequence, the differential conductance drops with the bias voltage, 
leading to the zero-bias peak, see the dashed line in fig. [2(c). This is thus the nonequilibrium 
accumulation of spin S = 1 which is responsible for the maximum in G at low bias voltage. 
Because the dot is coupled symmetrically to the leads, there is no spin accumulation in the 
parallel configuration and the differential conductance exhibits a smooth parabolic dependence 
on the bias voltage, see the solid line in fig. Etc). 

In the case of Agx > 0, at equilibrium the dot is occupied by a triplet. Therefore, at 
low bias voltage, \eV\ < |Agx|, transport is mediated by S = I states. When the voltage is 
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increased, at \eV\ ~ | Agx|, the occupation of the singlet state |0, 0)2 increases and all the four 
states start to participate in transport. This leads a step in differential conductance which is 
clearly visible in both magnetic configurations, see fig.|2d). The width of the transport region 
where cotunneling through the singlet state is suppressed is given by 2|Ast|- As one can see 
in the figure, in the parallel configuration there is one minimum, while in the antiparallel 
configuration there are two minima separated by a zero-bias peak. The mechanism leading 
to the peak at zero bias in the antiparallel configuration is the same as in the aforementioned 
case of Ast = 0. Spin accumulation induced in the dot, see fig. |3| leads to a decrease in 
differential conductance when increasing the bias voltage from V = 0. 

The nonequilibrium spin accumulation in the antiparallel configuration as a function of 
the bias voltage is displayed in fig. |3| In the case of Ast > 0, spin accumulation is always 
present for V 7^ 0, while for Ast < 0, spin accumulation appears when \eV\ « |Ast|, i.e. at 
the onset of cotunneling through triplet states. As pointed previously, spin accumulation is 
responsible for the maximum in differential conductance. Thus, one could expect that some 
maxima in G should also appear in the case of Ast < 0, where spin accumulation is induced 
for \eV\ > |Ast|, see fig. EH This is indeed visible in fig.[2a) and (b) where two small maxima 
develop at \eV\ « |Ast|> symmetrically with respect to the zero bias. 
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Fig. 4 - (color online) Tunnel magnetoresistance as a function of the bias voltage for Ast = 
—5, —2.5, 0, 2.5F, as indicated in the figure. The other parameters are the same as in fig.|21 

Another feature visible in fig. [5] is the difference between conductance in the parallel and 
antiparallel magnetic configuration. Conductance in the parallel configuration is generally 
larger than that in the antiparallel configuration. This is due to the asymmetry between 
the couplings of the dot levels to the spin majority and spin minority electron bands in the 
antiparallel configuration [7]. This difference gives rise to tunnel magnetoresistance, TMR = 
(Ip — IapV^ap, where ip (Zap) is the current flowing through the system in the (anti)parallel 
configuration [7,9]. The bias voltage dependence of TMR is shown in fig. 0] As one can 
see, TMR displays a nontrivial dependence on the ground state of the dot. When the dot is 
occupied by a singlet, Ast < 0, TMR at low bias displays a maximum plateau. With the 
same assumptions as made when deriving eq. (5), one finds that the maximum in TMR for 
A ST < and \eV\ < |A ST | can be approximated by, TMR S=0 = 2p 2 /(l - p 2 ). This value 
corresponds to the Julliere TMR, which is characteristics of single tunnel junction [9]. In our 
case, it results from the fact that in this transport regime current flows due to non-spin-flip 
cotunneling. On the other hand, for Ast > 0j the TMR exhibits a minimum at zero bias, 
which can be seen in fig.^for Ast = and A S t = 2.5T. 
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Summary. We have considered spin-polarized transport through a doubly occupied 
two-level quantum dot coupled to ferromagnetic leads in the cotunneling regime. In this 
case, due to finite exchange interaction, cotunneling current is mediated by singlet and triplet 
states of the dot. We have shown that transport strongly depends on the ground state of the 
system. If the dot is occupied by a singlet at equilibrium, G exhibits a minimum at low bias, 
whose width is given by the strength of splitting between the singlet and triplet states. In this 
transport regime the TMR is given by the Julliere value, indicating that transport is coherent. 
If triplet is the ground state, a zero-bias peak evolves in differential conductance for antiparallel 
configuration. The zero-bias maximum results from the nonequilibrium accumulation of spin 
S = 1 in the dot. The maximum is accompanied by a minimum in TMR. The different 
behavior of differential conductance and TMR may help in determining the ground state of 
the dot as well as the splitting between the singlet and triplet states in the dot. 

* * * 
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